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Abstract 



We describe a perturbation expansion for the energy and wave function of a weakly 
bound particle in a short-range potential in one space dimension. 
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1 Introduction 



In this article, we discuss a method for calculating energies and wave functions for weakly bound 
particles in short range potentials in one dimensional nonrelativistic quantum mechanics. We usually 
think of the formation of a bound state as a non-perturbative process, however, we will see that for a 
weakly bound state, with a wave function that extends far outside the region in which the potential 
is important, we can calculate bound state properties in perturbation theory. In one dimension, 
the perturbation expansion completely determines the energy and wave function of the state. The 
expansion is quite simple (closely related to the Born series) and it may well appear somewhere in 
the literature, but we have been unable to find it. 

A particle moving in a generic short range potential has two natural length scales: a, the range 
of the potential, and the size of the state, related to the binding energy, Eq, of the system 
by k = y/— 2mE ^ In turn, these scales define a dimensionless parameter, Ka. A weakly bound 
state, one for which Ka < 1, inherently has a small dimensionless parameter which could serve as 
the basis for a perturbative description of the wave function. In the following, we show that in 
one dimension we can construct a power series in na for both the energy and (in a sense) the wave 
function, assuming that na <ti l.@ 

The perturbative expansion that emerges is a series of n coupled integrals with n powers of 
the potential and n — 1 coupling terms at each order of the power series. These terms manifestly 
respect translation invariance — if the potential is translated in space, the wave function moves 
correspondingly. For simple potentials, the analytic behavior of the expansion for the bound states 
can easily be extracted to several orders. We illustrate the procedure for a classic square well 
potential in order to demonstrate the agreement of the perturbative approach with the known form 
of the ground state. Less simple potentials can be handled numerically. 

Our goal is eventually to extend this perturbative approach to three dimensions and from there to 
fully relativistic systems. The extension of these one dimensional results to three is not as straight- 
forward as might be hoped. The particle's wave function varies significantly where the potential is 
large, making the unraveling of the leading na behavior considerably more difficult. Moreover, the 

throughout this article, we use units in which % = 1. 

2 The situation in three dimensions is much more complicated. This condition is replaced by an integral condition 
on the potential — the small parameter then is the deviation of the actual potential from this condition. 
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form of the free particle Green's function in three dimensions, unlike its one dimensional analogue, 
is singular in its first term when expanded in powers of na. In this note, however, we set aside these 
difficulties and turn our attention to the particle in one dimension, where the expansion is simple 
and elegant. 



2 The Perturbative Expansion 

In the limit that the potential is very short range, we expect that it can be approximated by a 
5-function: 

V{x)^v5{x) (2.1) 

where 

v = J dxV(x). (2.2) 
The 5-function potential has a single bound-state, of the form 

ip(x) oc e- K|x| (2.3) 

where 

k = —mv (2-4) 

We seek a systematic expansion in which ( j2.1|) -( |2"/ID is the first term. Note that because k oc V(x), 
we expect that an expansion in powers of Ka will be an expansion in powers of V as well. A simple 
expansion of the Fourier transform of V(x) in powers of the momentum is not suitable because the 
higher order terms, derivatives of 5-functions, are too singular. Instead, we will exploit the relation 
between the ground state of the 5-function potential and the free particle Green's function, G K : 

G K (x) = e -^-. (2.5) 
In integral form, the Schrodinger equation reads 

ip(x) = - J dx' G K (x - x') 2mV(x')i)(x') (2.6) 
where k is related to the energy of the state, E , by 



-2mE . (2.7) 
The weakly bound regime occurs when na 1 where a is the range of the potential. 



If the argument of the Green's function is small compared to 1/k, we can legitimately expand 
the exponential in equation [2~5i In general, in (|2.6|) , this is not true. However, if we multiply ( |2.6|) 
by V(x), then we get an integral equation for the product, V(x)ip(x)\ 

V(x)tfj(x) = -V(x) J dx' G K (x - x') 2mV(x')ip(x > ) . (2.8) 

In ( p.8|) , the short range of the potential guarantees that can do this expansion because the x — x' 
in the Green's function is sandwiched between V(x) and V(x'), so \x — x'\ cannot get much larger 
than a. 

By adding and subtracting 1 from the exponential in G K , we can rewrite ( [2.8|) as 



V(x)iJj(x) 



Jdx'V{x')^{x') 



(2.9) 



mV(x) 



J dx' 



-k\x— x'\ 



V{x')^{x'). 



Since the principal contribution to the second integral comes only from the region x, x' < a, the 
second term on the right hand side of fl2.9p is explicitly suppressed by a factor of na relative to 
the first. The expression / dx' V(x')il){x') appears as a factor in all the terms when V(x)ip(x) is 
expanded to an arbitrary order so it is to convenient to work from the start instead with the 
normalized function 

V(x)ip(x) 



X{x) 



satisfying 

In terms of x ? we have 



Jdx'V(x')i{;{x') 
dxx(x) = 1 . 



mV(x) mV(x) 



dx' 



- 1 



K K 

This expression for x(x) serves as the basis for our perturbation series. 



(2.10) 
(2.11) 

(2.12) 



2.1 The Energy 



From the power series for x( x ) m equation |2.12| , we can obtain an expression for the bound state 
energy by applying the normalization condition, ( |2.11| ). For example, simply neglecting the second 
term in ( p.l2j ) and imposing ( |2.11|) yields 



/ 



dx 



mV(x) 



(2.13) 



which is equivalent to ( |2.4| ). Indeed, for a ^-function potential, the second term in ( [2.12| ) vanishes, 
thus ( p,13|) is exact, and we reproduce the solution (|2.3|) . 

Iterating (|2~T2l) once and neglecting terms of order m 3 V 3 yields 



mV(x) 



2mV(x) / dxi 



2k 



mV(xi) 



The leading contribution to \ is 



mV(x) m 2 V(x) 



dxi \x — x\\ V{x\) + 



K K 

for which the normalization condition determines the energy to be 

n = —m / dxV(x) — m 2 / dx\dx2\x\ — a^l V(xi)V(x2) + 



(2.14) 



(2.15) 



(2.16) 



It is straightforward to continue the expansion to higher order. The expression up to fourth order 
for k is given in the appendix, in equation |A.24j . 

Although the terms begin to proliferate as we expand to higher orders, the set of integrals that 
contribute with an overall coefficient of m n can be described. The set consists of all inequivalent 
integrals that contain n integrations, n factors of the potential, and n — 1 factors of the form 



x i+i| mj where the sum of the exponents satisfies 2j=i f^j 
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1. Since the \xj — 



arises from our expansion of the Green's function, 
general n-th order term looks like 



it is accompanied by a (l/rrij\). Thus, 



- m n s 



II dx i V( 



\i=l 




(2.17) 



There only remains the factor s to explain. In demanding that the integrals be inequivalent, we 
have slightly undercounted, and the s corrects for this omission. A specific contribution to (|2 . 1 7|) 
often contains several connected integrals and with each such piece written as 



dx\ • • • dxi~< n \x\ — X2 \ mi \x2 — X3 



I l»2 



\x k _ 1 -x k \ m ^V(x 1 )---V(x k ). 



(2.18) 



s — 1, then, counts the number of non-trivial ways of rearranging the rn^'s within this connected 
factor which yield different integrands but the same integral. Similarly, if our term ( 2.1 7|) contains 
several connected integrals, the total s is the product of these symmetry factors for each connected 
component. For example, at fourth order, shown in appendix equation ( |A.24 ), we have a term of 
the form, 

dx V(x) / dx\ dx2 dx3 \x% — X2\ 2 \x2 — Xs\ V (xi)V (x2)V (x 3 ) . (2-19) 



The exponents of the second connected piece could be switched so that the integrand has \x\ 



%2 1 \x2 — x 3 1 instead of | X\ — X2 I \x2 — X3 |; but both integrals are equivalent. Thus for (|2.19|) , s = 2 
(which cancels the 1/2! from the Green's function expansion). 



2.2 The Wave Function 

The expression for x( x ) to any order contains an explicit factor of V(x). One might think, from 
( |2.10| ), that we could simply cancel this factor and get the wave function up to a normalization 
factor. This is not quite right, because the resulting series only converges for small x. However, we 
can construct the wave function for all x using fl2.6| ), which can be written as 

ifj{x) = N J dx'e- Klx ~ x ' 1 x(x') (2.20) 

where N is a normalization factor. 



3 Example: The Square Well Potential 



For the square well potential, 



V(x) = < 



—Vq for \x\ < a 



for |x| > a 
our expansion is in powers of the dimensionless quantity, 

2mVoa 2 . 



If we apply ( A. 24 ), we find 



kcl = (2mV a 2 ) - l(2mV a 2 ) 2 + l(2mV a 2 
-f|(2m\/ a 2 ) 4 + O ({2mV a 2 ) 5 ) + ■■■ . 



(3.21) 



(3.22) 



(3.23) 



which agrees with what we get from a Taylor expansion of the exact value of n. 

Rather than show the rather complicated formulae for the expansion of the wave function, we 
will display, in figures 1-4, plots of the exact wave function along with the first four approximates, 
for various values of 2mVoa 2 . Evidently, for 2mVoa 2 = 0.3, shown in figure 1, the expansion for the 
wave function converges rapidly, while for 2mVoa 2 = 0.9, shown in figure 4, the expansion is not 
working at all. 



6 



Acknowledgements 

We are very grateful to David Kaplan for many interesting conversations. 

A Appendix: Expansion Formulae 

This appendix lists the energy and x expansions to several orders in mV; the structure of the 
collection of terms present at a particular order in the na expansion has been described in section 
2. 



k = —m J dx V(x) — m 2 J dx\ dx 2 \x\ — x 2 \ V(xi)V(x 2 ) 



-m 



J dx 1 dx 2 dx 3 \x 1 - x 2 \\x 2 - x 3 \ V \xi)V \x 2 )V \x 3 ) 



— m 



+h[Jdx V(x) )\Jdx 1 dx 2 \ Xl - x 2 \ 2 V (x^V \x 2 ) 



J dx\ dx 2 dx 3 dx4 \xi — x 2 \ \x 2 — x 3 \ \x 3 — x 4 | V{x\)V{x 2 )V{x 3 )V{xi) 



+ dxV(x) S j ^fdxi dx 2 dx 3 \xi - x 2 \ 2 \x 2 - x 3 \ V(x 1 )V(x 2 )V(x 3 ) S j 
+l^J dxi dx 2 \xi - x 2 \ 2 V(xi) V(x 2 )j dxi dx 2 \x 1 - x 2 \ V(x 1 ) V(x 2 ) 
+|( / dxV(x) I I / dx 1 dx 2 \x! - x 2 \ 3 V '{xi)V \x 2 ) 



(A.24) 



+0(m 5 V 5 ) 
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—kx(xi) = mV(xi) + m 2 J dx 2 \x\ — x 2 \ V(x\)V (x 2 ) 



+m 6 



f dx 2 dx 3 \x 1 - x 2 \\x 2 - x 3 \ V \xi)V (x 2 )V (x 3 ) 



+m 



+\{jdx 2 \x 1 -x 2 \ 2 V{x l )V{x 2 )^ ^JdxV(x)^ 



J dx 2 dx 3 dx 4 \xi - x 2 \\x 2 - x 3 \\x 3 - x±\ V(xi) V(x 2 ) V(x 3 ) V(x 4 ) 



+\(yS dx 2 dx 3 \x x -x 2 \\x 2 - x 3 \ 2 V(x 1 )V(x 2 )V(x 3 ) S j (^JdxV(x) 
+l^J dx 2 dx 3 \xi - x 2 \ 2 \x 2 — x 3 \ V(x 1 )V(x 2 )V(x 3 ) S j dxV(x) 
+ \ ( / dx 2 \ Xl - x 2 \ 2 V(x 1 )V(x 2 ^j dx 2 \ Xl - x 2 \ V{x 1 )V{x 2 ) S j 
+\{^dx 2 \x 1 -x 2 \ z V{x 1 )V{x 2 )^ ^JdxV(x)^ 



+0(m 5 V 5 ) 
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Figure 1: 2mV a 2 = 0.3 
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Figure 2: 2mV a 2 = 0.5 
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